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The viscoelasticity of a two-dimensional liquid strongly-coupled dusty plasma is studied exper- 
imentally, without macroscopic shear. Positions and velocities of the dust particles, measured by 
video microscopy, are used as the inputs to the generalized Green-Kubo relation to obtain the com- 
plex viscosity rj(uj). The real part of (which corresponds to dissipation) diminishes gradually 
with frequency, while the imaginary part (which corresponds to elasticity) is peaked at a frequency 
below the 2D dusty plasma frequency. The viscoelastic approximation is found to accurately de- 
scribe the 2D experimental results for f](ui), yielding the Maxwell relaxation time tm = 0.10 s. 
Results for r?(a>) are compared to 2D molecular dynamics Yukawa simulations and to a previous 
experiment that was performed using an oscillating macroscopic shear. 
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I. INTRODUCTION 

Liquids behave as a viscous continuum in the hydro- 
dynamic regime of large distances and long times, but in 
the kinetic regime where distances and times are compa- 
rable to molecular scales, liquids can also exhibit clastic 
effects Q. Viscous effects correspond to energy dissipa- 
tion, while clastic effects correspond to energy storage. 
When they occur at the same time, as they can in a 
liquid in the kinetic regime, the phenomenon is called 
viscoelasticity 0. 

To characterize the viscoelasticity of a liquid quantita- 
tively, one can generalize the concept of viscosity, so that 
it is dependent on both frequency cj and wave number 
k @, 0|- It is common, however, to report a simplified 
quantity that depends only on frequency rj(u>), which is 
a complex function: 



T](ui) = T)'(w) - iT)"(u). 



(1) 



The real part rf(ui) corresponds to the liquid-like prop- 
erty of viscous dissipation, while the imaginary part 
t]"(uj) corresponds to the solid- like property of elastic- 
ity [5kZ|. By determining r)(to), one can compare the 
effects of energy storage and dissipation of a liquid for 
a particular time scale A simple model of the 

frequency-dependent viscosity is predicted by the vis- 
coelastic approximation [l[ which has two parameters: 
the Maxwell relaxation time tm and the static viscosity 
rjo, as will be discussed in Sec. II A. The latter quantity is 
the same as t){uj) in the hydrodynamic limit of ui — > 
Aside from t](oj), one can also characterize viscoelasticity 
using the wave- number-dependent viscosity r](k), as has 
been done in molecular dynamics (MD) simulations of 
liquids consisting of particles that interact through vari- 
ous potentials [7Hllj and also in a two-dimensional (2D) 



dusty plasma experiment [11 



In this paper, we will investigate viscoelasticity of 2D 
liquids by performing experiments with video microscopy 
observations of particle motion in a liquid. There are 



several 2D physical systems that allow this kind of direct 
observation of individual particle motion. These include 

granular materials |1 4[| , a 



colloidal suspensions [12 



13] 



Wigncr lattice of electrons on a liquid helium surface [15| , 
ions confined magnetically in a Penning trap [lif , vortex 
arrays in the mixed state of type-II superconductors [13] , 
and strongly-coupled dusty plasmas levitated in a single 
layer [HEl- 

A dusty plasma UH)], sometimes termed complex 
plasma, is a partially ionized gas containing solid micron- 
size particles, which we will refer to in this paper as 
dust particles. Due to their large length scales and 
slow time scales [2lJ, dusty plasmas allow imaging by 
video microscopy and tracking of individual particle mo- 
tion (2^1 ■ Dust particles are negatively charged, with typ- 
ically thousands of elementary charges. As a result, an 
ensemble of mutually repulsive dust particles is strongly 
coupled 2(| 22j. Dusty plasma experiments can be per- 



formed with all the dust particles located in a single hor- 
izontal layer that is perpendicular to the ion flow; under 
these conditions, the collection of dust particles behaves 
essentially as a 2D system, with interparticlc interactions 
that can be modeled through a repulsive Yukawa poten- 
tial 2H 2^ with a screening length Xrj. Using various 
heating methods [soI-EmI]. the kinetic temperature [HI 
can be adjusted so that the suspension of dust parti- 
cles behaves like a liquid or solid. In previous dusty 
plasma experiments, elasticity in solids |35( and viscosity 
in liquids |36|-|39j have been studied by applying shearing 
stress using laser beams. 

We will obtain rjiuS) for a 2D dusty plasma using exper- 
imental data in the generalized Green-Kubo relation, as 
discussed in Sec. II C. In our literature search, we have 
not found any previous report, for any substance, of a 
similar use of the generalized Green-Kubo relation with 
an input of experimental data. The experimental data 
we use is from a 2D dusty plasma experiment that was 
performed earlier in our lab [4ol ]. to maintain a 2D liq- 
uid under well-controlled steady conditions. These condi- 
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tions are desirable for using the generalized Green-Kubo 
relation. 

We also test the expressions of the viscoelastic approx- 
imation, which is a standard theory for simple liquids 
that has previously been tested experimentally for liquid 
argon and other 3D liquids but not, to the best of 
our knowledge, for 2D liquids or dusty plasmas. We per- 
form this test by fitting the equations for the viscoelastic 
approximation to the rj(u)) measured in two experiments, 
from our lab [4(J and the lab of Hartmann et al. fl8j . 
both using 2D dusty plasmas. 



In the viscoelastic approximation [![, rj(uj) can be ex- 
pressed as 



+ 1/tm ' 



(2) 



where G^ is an instantaneous modulus of rigidity, and 
Tm is the Maxwell relaxation time, tm = Vo/Goo- The 
Maxwell relaxation time tm characterizes the transition 
between viscous behavior at low frequencies {ujtm "C 1) 
and elastic behavior at high frequencies (ujtm 3> 1) 01- 
The real and imaginary parts of Eq. ([2]) are: 



II. FREQUENCY-DEPENDENT VISCOSITY 

In this section, we review the principles and methods 
for the viscoelastic approximation, and methods of ob- 
taining the frequency-dependent viscosity rj(oj). We will 
use these principles and methods in Sec. IV. 



A. Viscoelastic approximation 

The frequency-dependent viscosity r](uj) characterizes 
how both viscous and elastic effects vary with time scales 
or frequencies. The two parts, rj'(u)) and 7]"(u), corre- 
spond to viscous dissipation and elasticity, respectively. 
As tu — > 0, it is expected that the real part extrapolates to 
the usual static viscosity 770 , i.e., zero frequency viscosity. 




frequency co (arb. units) 

FIG. 1: (Color online). Sketch of the viscous and elastic 
parts of the frequency-dependent viscosity, Eq. (jTJ) , as in 
At low frequencies, motion is dominated by dissipation, cor- 
responding to the real part 77', while at high frequencies, it is 
dominated by elasticity, corresponding to the imaginary part 
V"- 



for the viscous part, and 



1 + T> 2 



(3) 



(4) 



for the elastic part. Equations (O and Q have previously 
been used to fit rj(uj) from a simulation of a 3D Yukawa 
liquid [41 1 . 

In Sec. IV, we will fit Eqs. ([3]) and ((4]) to our experi- 
mental measurements of rj(u>), with one free parameter, 
tm- We will determine tm for our experiment of [ioj ]. 
and compare it to values we determined from data for 
another experiment [l8| and our MD simulations. 



B. Rheometry to measure 77(0;) 

A rheometer is an instrument that imposes a macro- 
scopic shear to measure viscosity. A boundary contact- 
ing the fluid is rotated, causing a flow with a velocity 
gradient. The resulting shearing stress can be measured 
to determine the viscosity 77, using its hydrodynamical 
definition [ijj ]. This can be done either with a steady 
or oscillating rotation, to find 770 or t](oj), respectively. 
The imaginary part of rj(uj) is determined by measuring 
a phase shift (i.e., delay) between the applied modulation 
and the resulting shearing stress. 

Dusty plasmas are not suited for use with rheometers, 
however, because the dust particles do not contact any 
solid surface. Therefore, experimenters must determine 
viscosity some other wa y, s uch as applying a shearing 
stress using laser beams 3, 37 1, or by using the Green- 
Kubo relation [4(J as we will do in this paper. 



Figure 1 shows a sketch of a typical complex viscosity 
r](u)) for a liquid. Both the real and imaginary parts of 
r](ui) have been normalized by the static viscosity 770. As 
the frequency u increases from 0, the real part 7/ (w) de- 
cays from unity monotonically toward 0. The imaginary 
part 77" (w) increases from to its maximum and then 
decays gradually. 



C. Generalized Green-Kubo relation to measure t}{uj) 

The generalized Green-Kubo relation can be used to 
obtain the frequency-dependent viscosity. Previously it 
was used in simulations 3] • It can be used for conditions 
without macroscopic shear. Its inputs are the positions, 
velocities and potential energies of individual particles 
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as they move in their thermal motion. The frequency- 
dependent viscosity r](uj) is calculated in three steps. 
First, the shearing stress P xy (t) is calculated as 



N 



Pxy{t) = 



1 ^ Xijyij d$(rij) 



' y 



dr 



i j 



(5) 



Here, i and j indicate different particles, N is the total 
number of particles of mass m, r, = (x iy y^) is the position 



of particle i, r.j 



and &(rij) is the intcrparticlc potential energy. 

Second, the shearing stress autocorrelation function 
(SACF) is calculated as 



C v (t) = (P xy (t)P xy (0)). 



(6) 



For the experiment reported here, despite the efforts that 
were made, there was still a small nonzero average veloc- 
ity in the analyzed region, which results in a nonzero 
time-averaged P xy . To solve this problem, as in [40j . we 
subtract the time-average value from P xy (t), so that we 
use only the fluctuating portion to calculate the SACF. 

Third, the frequency-dependent viscosity t)(uj) can be 
calculated from the Laplace-Fourier transformation of 
C v (t), i.e., the generalized Green-Kubo relation: 



rj(uj) 



1 



Ak B T 



C v (t)e iuit dt, 



(7) 



where A is the area of the 2D system and T is its tem- 
perature. Equatio n ([71 ) can be derived by combining Eqs. 
(20) and (21) of |43( for a zero wavenumber, and it has 
been used previously for a 3D Yukawa simulation Q • The 
units of r\ are kgs -1 for a 2D liquid (unlike the viscos- 
ity of a 3D liquid, which has units of kgm _1 s _1 ). By 
specifying u = 0, Eq. (J7J is the same as the usual Green- 
Kubo relation for the static viscosity rjo (HI- Unlike the 
shearing stress in Eq. ([5]), 7](u>) in Eq. ([7j is independent 
of the system size. 

The required inputs to calculate P xy (t) are time se- 
ries of positions, velocities, and potential energies of 
each particle. For experiments that rely on video mi- 
croscopy, however, only positions and velocities are usu- 
ally recorded, not $. For the experiment reported here, 
the potential energy cannot be measured directly, so 
we calculate it from other measurements by assuming 
a model potential. We use the Yukawa potential, which 
has been shown experimentally [28| and theoretically [1^] 
to describe the conditions in a plane perpendicular to the 
ion flow at a distance of about Ad , as is the case for the 
single layer of dust particles in our experiment [40| . 

To obtain the frequency-dependent viscosity rj(io) us- 
ing Eq. ([7J with an input of experimental data, it is nec- 
essary to replace the infinity upper limit with a finite 
time. In Sec. IV, we will replace the upper limit of the 
integral in Eq. (J7J with the time that C n (t) crosses zero 



for at least two consecutive data points in the time series, 
as in [40| . This choice for the upper limit determines the 
smallest value of w at which we can measure r)(cd). As in 
Eqs. ([3]) and ((4]), we normalize r)(ui) by the static viscos- 
ity Vo, which we calculate using Eq. ([7J with uj = 0, as 
in libl. 



D. Review of r\{(J) in dusty plasmas 

The viscoclasticity of dusty plasmas has been the sub- 
ject of several theoretical studies 44-49j. It has also 
been studied using 3D Yukawa MD simulations 0] ■ Early 



experimental studies [50j, l5l| were followed by measure- 
ments of the wave-number-dependent viscosi ty r? (fc) [ll[ 
and the frequency-dependent viscosity r\{iS) [l8j. In the 
latter experiment [18j . a macroscopic oscillatory shear 
was applied to a 2D dusty plasma using a time-modulated 
laser manipulation. 

The generalized Green-Kubo relation Eq. (JTJ) that we 
use in this paper is based on the fluctuating shearing 
stress associated with the thermal motion of the dust 
particles, without any macroscopic shear. Equation ([7]) 
has been tested using 3D simulations Q where it was 
shown that it yields t\(uj) consistent with results from a 
simulation using a different method, with an oscillatory 
shear. Our use of Eq. (J7J is different because we use 
experimental data, and we do this for a system that is 
2D instead of 3D. 

In this paper, we use the generalized Green-Kubo re- 
lation with an input of experimental data to obtain the 
complex viscosity r\{uS). The experiment, which was re- 
ported in [io[, is reviewed briefly in Sec. III. Our calcu- 
lation differs from that of [4(J because the Green-Kubo 
relation that we use, Eq. ([7]), is the generalized form, 
so that we obtain the frequency-dependent viscosity in- 
stead of the static viscosity t\q. In Sec. IV we present 
our results for r\(iS) and tm, and we compare to the pre- 
vious experimental results of Hartmann et al. 18[, who 



used an applied oscillatory shear. In Sec. V we compare 
our experimental results to 2D MD simulations, with and 
without the effects of gas friction. 



III. EXPERIMENT 

In this section, we briefly review the experiment of (40| , 
which provides the data we will use in Sec. IV. Dust parti- 
cles, which were 8.1 /itm diameter polymer microspheres, 
were electrically levitated as a 2D horizontal layer in an 
argon rf plasma. Using a top-view camera operated at 
250 frame/s, we tracked w 2100 dust particles, yielding 



time series of their positions and velocities [52|, |53j . Four 



runs, each lasting 20.2 s, were performed with laser heat- 
ing so that the dust particles behaved collectively as a liq- 
uid. In this liquid, the dust particle kinetic temperature 
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was nearly uniform, with a value of T = 2.5 x 10 4 K. The 
"kinetic temperature" describes the motion of a dust par- 
ticle, not the temperature of the polymer material within 
the dust particle. We calculate the kinetic temperature 
from the mean-square velocity fluctuation as in [33 1. 

The layer of dust particles had the following param- 
eters. The interparticle distance, characterized by the 
Wigner-Seitz radius HH a = 0.35 mm, corresponded to 
an areal mass density of p = 1.1 x 10~ 6 kg/m 2 . A moving 
dust particle experienced a drag force due to the neutral 
gas characterized by a friction coefficient v$ = 2.4 s _1 . 
The charge on a dust particle was QJe= —6000, the 2D 
dust plasma frequency was tu p( i — 30 s" 1 , the cou- 
pling parameter was T — (Q 2 / Aire^a) / (UbT) = 68, and 
the screening parameter was n = a/Azj = 0.5. The values 
of Q and n have uncertainties estimated to be ±10 % and 
±20 %, respectively. Further details of this experiment 
can be found in [401 ] . 

The result for the shearing stress autocorrelation func- 
tion (SACF) for one experimental run is shown in Fig. 2. 
We calculate the shearing stress autocorrelation func- 
tion (SACF) using Eq. ©, after obtaining P xy (t) using 
Eq. ([5]). The SACF generally has a trend of decay, start- 
ing with its maximum value at t = 0. The frequency 
content within this decay is important in the determina- 
tion of T)(ui). 
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FIG. 2: (Color online). The shearing stress autocorrelation 
function (SACF), i.e., C v , calculated using an input of experi- 
mental data from one experimental run. The vertical axis, the 



SACF, is normalized by (AksT [ pa 2 uj^ d ) 



while the horizontal 



axis, time, is normalized by w. 



pd 



In Sec. IV we will analyze data from four runs to ob- 
tain the frequency dependent viscosity rj(uj) using Eq. ([?]). 
This result for ?/(w) will have an uncertainty of about 13% 
due to the uncertainties in Q and k. 



IV. RESULTS 
A. Frequency-dependent viscosity t](lj) 

In Fig. 3 we present our chief result: the frequency 
dependence of the complex viscosity, 77(01). These data 
points are obtained using the generalized Green-Kubo 
relation, Eq. (JT)), after the SACF is calculated, for each 
run. We also obtain the zero-frequency static viscosity 
770 for that run by using Eq. ((TJ) with ui = 0. Repeating 
for all four runs, we obtain the data for r](uj) shown in 
Fig. 3. 
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FIG. 3: (Color online). The frequency-dependent viscosity 
rj(to) obtained using an input of experimental data in Eq. (0. 
These results are shown as data points for four experimental 
runs. The smooth curves are the result of a simultaneous fit 
using the viscoelastic approximation, Eqs. and ([J). This 
fit yields the Maxwell relaxation time, tm = 0.10 s = 3.0u>~ d . 
The vertical axis is normalized by the static viscosity 770 for 
each run. 



The trends seen in our results for t]'(ui) and 77" (w) in 
Fig. 3 for our 2D liquid are generally the same as for 3D 
liquids. The real part T)'(uj) (which corresponds to dissi- 
pation) decays with increasing frequency. The imaginary 
part 77" (w) starts at zero for ui = 0, then has a peak 
followed by a decay at large lj. 

To allow a comparison to other experiments and simu- 
lations, our results for the frequency-dependent viscosity 
77(0;) in Fig. 3 are presented in normalized units. Vis- 
cosity is normalized by the static viscosity 770, which was 
determined separately for each run. Frequency is nor- 
malized as u/ujpd, where uj p d is the 2D dust plasma fre- 
quency. The scatter in our results for rj(uj) in Fig. 3 arises 
from noise in the SACF, which is believed to be statis- 
tical noise due to the finite amount of data used in a 
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calculation for a single experimental run. 

B. Test of the viscoelastic approximation 

1. Without macroscopic shear 

We now test the viscoelastic approximation, using 
the results of the frequency-dependent viscosity obtained 
from the experiment of [40j], which had no macroscopic 
shear. We fit the expressions for the viscoelastic approx- 
imation, Eqs. and to our experimental results 
for the frequency-dependent viscosity data, yielding the 
smooth curves in Fig. 3. In performing this fit, we com- 
bined the data for all four runs in Fig. 3. This fit has 
a single free parameter, which is the Maxwell relaxation 
time tm- 

We find tm = 0.10 s, which in dimcnsionless units is 
3.0 The physical significance of the Maxwell relax- 
ation time is that elastic behavior dominates at shorter 
times t <C tm while viscous dissipation dominates at 
longer times t tm- 

2. With macroscopic shear 
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FIG. 4: (Color online). The frequency-dependent viscosity 
t](oj) from the experiment of Hartmann et al. which was 
performed using an oscillatory shear. Here, the viscosity has 
been normalized by pa?u) v a, differently from Figs. 3 and 5. 
The smooth curves are our fit using the viscoelastic approxi- 
mation, yielding tm = 3.1 w^ 1 and 770 = 0.16. Another time 
scale of interest corresponds to gas friction, which in this ex- 
periment is quantified by Vf = 0.04aj p[ j. 

For comparision, we also test the viscoelastic approx- 



imation, this time using data from the experiment of 
Hartmann et al. [l8j, which had a macroscopic oscillat- 
ing shear. That experiment was similar to ours in sev- 
eral ways. Both experiments used a single layer suspen- 
sion of melamine-formaldehyde microspheres immersed 
in a radio-frequency argon plasma, and the dust parti- 
cles were tracked using video microscopy. 

The two experiments used different laser manipulation 
schemes and different approaches to determine i](uj). We 
used a uniform laser heating and the Green-Kubo ap- 
proach, which requires that the dusty plasma have steady 
uniform conditions, while Hartmann et al. applied an os- 
cillatory shear at an adjustable frequency and obtained 
the viscosity as the ratio of the shearing stress P xy and 
the applied shear rate. The Green-Kubo approach tends 
to generate data most easily at high frequencies, because 
we replace the infinite time limit in Eq. (JT]) with a finite 
time, as mentioned in Sec. II. On the other hand, the 
data reported by Hartmann et al. were mainly for lower 
frequencies. Thus, the two methods are complementary 
to each other. 

Some of the experimental conditions were different. 
In (40| the dust-particle diameter was nearly twice as 
large as in [l8j], and the gas pressure was slightly higher 
in |40j . The screening length A d was approximately twice 
as large in even though the particle spacing was 

nearly the same, so that the screening parameters k were 
different, 0.5 for but 1.2 for [H]. 

We find that the viscoelastic approximation is capable 
of describing both experiments. This is seen in Fig. 4, 
where we see that the fit of the viscoelastic approximation 
Eqs. ([3]) and Q passes through the data points of Hart- 
mann et al. 1 1 81 ] . The values of the fit parameters tm and 
770 for the two experiments are not expected to match, 
since the experimental conditions were different. Never- 
theless, we find that both experiments had the same nor- 
malized static viscosity [!o| 770/ pa 2 Lu p d — 0.16. We also 
find that both experiments had nearly the same normal- 
ized relaxation time: tm = 3.0 to^} and tm = 3.1 uiZj, 
for the experiments of [40j and [18J, respectively. 



V. SIMULATIONS 

To demonstrate that the viscoelasticity that we ob- 
served for our dusty plasma experiment is the result of 
interactions of the dust particles, we performed MD simu- 
lations. These simulations used only simple physics, with 
an integration of the equations of motion of dust particles 
interacting through a binary Yukawa potential. Using 
the particle positions, velocities and potentials recorded 
in the simulations, we calculated r}'{uj) and t]"(uj). We 
then fit Eqs. ([3]) and (j4|) to these results to obtain tm- 

Two simulation methods were used: the friction- 
less MD and the Langevin MD methods, as described 
in 0] . Both used a binary interparticle interaction with a 
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FIG. 5: (Color online). Simulation results for a 2D Yukawa 
liquid. The frequency-dependent viscosity is accurately de- 
scribed by the viscoelastic approximation, shown by the 
smooth curves. The simulation used the same values of T 
and k as for our experiment of |4(j, and it was performed two 
ways, using the (a) frictionless and (b) Langevin methods. 



Yukawa pair potential for particles that were constrained 
to move in a 2D plane. For both simulations, we used 
N = 4096 particles in a rectangular box with periodic 
boundary conditions. Trajectories r;(i) were found by in- 
tegrating equations of motion for all particles in MD sim- 
ulations [H, with the integration time step of 0.019 c^Q 1 - 
The simulation data trajectories were recorded for a du- 
ration 68 000 after a steady state was reached. This 
duration was more than 100 x longer than one experimen- 
tal run, so that the simulation has better statistics. The 



dimensionless parameters T and k were chosen to match 
the experimental values of T — 68 and k = 0.5 llj. In 
our Langevin MD simulation, we chose the gas friction to 
be the same as in the experiment, Vf = 0.08 Lo p d- Ewald 



summation [55j was not required because our 2D sim- 
ulation box is large enough compared to the screening 
length [Hi]. 

Simulation results for rf(w) and t]"{uj) are shown in 
Fig. 5. They exhibit the same trends as in the experimen- 
tal results. The Maxwell relaxation time, tm-, is found to 
be 3.4 o;~ d for the frictionless MD simulation, and 3.8 
for the Langevin MD simulation. These two values are 
both close to the experimental value tm = 3.0 for 
the same T and n. This agreement between the exper- 
iment and these simulations demonstrates that the vis- 
coelasticity observed in the experiment is mainly due to 
binary interactions among dust particles. 



VI. SUMMARY 

We experimentally determined the frequency- 
dependent viscosity rj(ui) for a 2D dusty plasma, using 
the generalized Green-Kubo relation. This measurement 
relies on the thermal motion of individual dust particles 
for experimental conditions without any macroscopic 
shear. We found that the real and imaginary parts of 
f](u>) show agreement with the viscoelastic approxima- 
tion, so that we were also able to determine the value of 
the Maxwell relaxation time, tm = 0.10 s = 3.0 oj~J for 
our experimental conditions. 

Wc performed 2D MD simulations, using the same pa- 
rameters as in the experiment, and wc found that the 
simulations predict rj(ui) and tm that agree with the ex- 
periment. This comparison to simulation indicates that 
the viscoelastic effects in the 2D dusty plasma, as quan- 
tified by rj(u)), are mainly due to the binary intcrparticle 
interactions among the dust particles. 

The generalized Green-Kubo relation that we used for 
conditions without externally applied shear yielded data 
for r](w) mainly for higher frequencies. In contrast, Hart- 
mann et al. [18| . who used an applied oscillatory shear, 
reported results for lower frequencies. We find that the 
viscoelastic approximation is capable of describing both 
experiments. 

This work was supported by NSF and NASA. 
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